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ABSTRACT 

Context. While there are a lot of investigations on the optical variability of flat-spectrum radio quasars (FSRQs), not much work has 
been done on the optical variability of steep-spectrum radio quasars (SSRQs). 

Aims. We investigate the optical variability in SSRQs. For comparison, the optical variability of FSRQs are also explored. 
Methods. We investigate the optical variability of 18 SSRQs and 15 FSRQs in the SDSS Stripe 82 region using SDSS DR7 released 
multi-epoch data covering about nine years. We determined the spectral index by fitting a powerlaw to SDSS ugriz photometric data, 
and explored the relationship between the spectral index and source brightness. 

Results. For all SSRQs studied, we detect variations in r band flux of overall amplitude between 0.22 mag and 0.92 mag in different 
sources. Eight of 18 SSRQs display a bluer-when-brighter (BWB) trend. In contrast, the variability amplitude of 15 FSRQs in r band 
ranges from 0.18 to 0.97 mag, and 11 of 15 sources show a BWB trend. We found an anti-correlation between the Eddington ratio and 
the variability amplitude in r band for SSRQs, which is similar to that in radio quiet AGNs. This implies that the thermal emission 
from accretion disk may be responsible for the variability in SSRQs, although the jet nonthermal emission cannot be excluded. In 
contrast, there is no correlation in FSRQs, which implies that the mechanisms of variability in FSRQs may be different from that in 
SSRQs. The optical continuum variability of radio loud broad absorption line quasars (BALQs) are investigated for the first time here 
on two sources with steep radio spectrum. Both radio loud BALQs show variations with amplitude of about 0.3 mag at r band. Their 
spectral variability all show a BWB trend. In a combined sample (18 SSRQs, and 44 FSRQs) of our sample with the FSRQs in Gu 
& Ai (201 1), we found a trend of a broader line width of broad Mg II emission line with steeper radio spectral index. It implies a 
disc-like broad line region (BLR) geometry may present in these quasars. In these 62 quasars, we found a ~ 57% source percentage 
showing BWB trend in FSRQs, whereas it is ~ 44% in SSRQs. 

Key words, galaxies: active - galaxies: quasars: general - galaxies: photometry 

1 . Introduction terpreted as a consequence of non-thermal emission from a rel- 

ativistic jet oriented close to the line of sight. In general, the 

Active galactic nuclei (AGNs) are characteristic of variability at variations in blazars are dorninate d by the jet emission. There 

almost all wavelength (e.g. Wnta 1996). Multiwavelength stud- ^ extensive explorations on the optical variability of blazars 

les of variations of AGNs have played important roles in ex- ( Ghisellini et al. 1997; Fan et al. 1998; Massaro et al. 1998 

plonng the physical conditions near the center of AGNs. Based Ghofjh et al 2000; clements & Carini 2001 . Raiteri et al 2001 

on the variability, the reverberation mapping method was de- vnlata et aL 2 002; Vagnetti et al. 2003; Wu et al 2005, 2007 

■ veloped to calculate the size of broad line region (BLR). From Gu et al 2006; Hu et al 2006; p oon et al 20Qg . Rani et al 

the photoiomzation model, the flux variations of BLR follows 2010; Gu & Ai 2QU) while it is generally accepted that the 

the flux variation of continuum ionization emission with certain nontherma l emission from relativistic jet oriented close to the 

time lag, which corresponds to the size of BLR (Peterson 1993). line of sight dominate the optical continuum, the situation seems 

Moreover, the disc-jet connections can be investigated by us- more complicated in FSRQs. It is not clear whether the nonther- 

ing the long term multiwavelength variability monitonngs (e.g. mal jet emission plays the main role in the optical variability. 

Chatterjee et al. 2009, 2011). Evidence of thermal emission in FSRQs has been detected in 

Intrinsically, the variations of AGNs are generally caused several cases usually during low activity states ( e .g. 3C279, Pian 

by the physical variations in jet and accretion disc. However, et al i 999; 3C273, Grandi & Palumbo 2004; 3C454.3, Raiteri 

the contribution of each component in the observed variability et al 2007) A recent work shows mat the optical variability in 

varies from source to source. Blazars, including BL Lac ob- the FSRQ PKS i 5 i -089 is mainly due to the thermal emission, 

jects and flat-spectrum radio quasars (FSRQs), are the most ex- and only during ma j or flares a contribution from the jet is seen 

treme class of active galactic nuclei (AGNs), characterized by (D' Ammando et al. 201 1). Moreover, the redder-when-brighter 

strong and rapid variability, high polarization, and apparent su- trend found in several FSRQs implies that the thermal emission 
perluminal motion. These extreme properties are generally in- 
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plays important role in the optical variability (e.g. Gu et al. 2006; 
Rani et al. 2010; Gu & Ai 201 1). 

Radio-quiet AGNs has also been explored in terms of the 
long-term or short-term variability (e.g. Stalin et al. 2004, 2005; 
Gupta & Joshi 2005; Ai et al. 2010). The accretion disc instabili- 
ties may likely explain the optical microvariability in radio quiet 
AGNs (e.g. Gopal-Krishna et al. 1995, see also Gopal-Krishna 
et al. 2003 for alternative scenario). In terms of the long-term 
variability, the change of accretion rate is used to explain the 
optical variations in radio quiet AGNs (e.g. Li & Cao 2008). 
The properties of steep spectrum radio quasars (SSRQs) are in- 
termediate between FSRQs and radio quiet quasars. SSRQs are 
usually lobe-dominated radio quasars, with radio lobe emission 
dominate over the radio core emission. Their jets are viewed at 
larger angles than blazars. Therefore, the beaming effects from 
jets should not be severe (see e.g. Liu et al. 2006), and the jet 
emission is not expected to dominate at optical bands (e.g. Gu 
& Ai 2011). However, the variations in SSRQs could still be 
a mixture of being caused by jet and accretion disc, since they 
are powerful radio emitter. The optical variations of SSRQs are 
largely unknown, and there are only few explorations on varia- 
tions in SSRQs (e.g. Stalin et al. 2004, 2005). While the color 
variations and/or spectral variations of blazars are extensively 
investigated (e.g. Ghisellini et al. 1997; Fan et al. 1998; Massaro 
et al. 1998; Ghosh et al. 2000; Clements & Carini 2001; Raiteri 
et al. 2001; Villata et al. 2002; Vagnetti et al. 2003; Wu et al 
2005, 2007; Gu et al. 2006; Rani et al. 2010; Gu & Ai 2011), 
it has rarely been done for SSRQs. In this work, we investigate 
the optical variability for a sample of SSRQs, as well as their 
spectral variations. 

The layout of this paper is as follows: in Section 2, we de- 
scribe the source sample; the variability results are outlined in 
Section 3; Section 4 includes the discussion; and in the last 
section, we draw our conclusions. The cosmological parame- 
ters H Q = 70 km s" 1 Mpc" 1 , fi m = 0.3, and Q A = 0.7 are 
used throughout the paper, and the spectral index a is defined 
as f v oc v~ a , where f v is the flux density at frequency v. 



2. Sample selection 

2.1. Quasars in Stripe 82 region 

Our initial quasar sample was selected as those quasars both in 
the SDSS DR7 quasar catalogue (Schneider et al. 2010) and 
Stripe 82 region. The SDSS DR7 quasar catalogue consists of 
105,783 spectroscopically confirmed quasars with luminosities 
brighter than M, = -22.0, with at least one emission line having 
a full width at half-maximum (FWHM) larger than 1000 km s _1 
and highly reliable redshifts. The sky coverage of the sample is 
about 9380 deg 2 and the redshifts range from 0.065 to 5.46. The 
five-band (u, g, r, ;', z) magnitudes have typical errors of about 
0.03 mag. The spectra cover the wavelength range from 3800 to 
9200 A with a resolution of 2000 (see Schneider et al. 2010 for 
details). The Stripe 82 region, i.e. right ascension a = 20 h - 4 h 
and declination 6 = -1°.25 - +1°.25, was repeatedly scanned 
during the SDSS-I phase (2000 - 2005) under generally photo- 
metric conditions, and the data are well calibrated (Lupton et al. 
2002). This region was also scanned repeatedly over the course 
of three 3-month campaigns in three successive years in 2005 - 
2007 known as the SDSS Supernova Survey (SN survey). The 
multi-epoch photometric observations therefore enable us to in- 
vestigate the optical variability of the selected quasars. 



2.2. Cross-correlation with radio catalogues 

In this paper, we define a quasar to be a SSRQ according to 
its radio spectral index. Therefore, we cross-correlate the initial 
quasar sample with the Faint Images of the Radio Sky at Twenty 
centimeters (FIRST) 1.4-GHz radio catalogue (Becker, White 
& Helfand 1995), the Green Bank 6-cm (GB6) survey at 4.85 
GHz radio catalogue (Gregory et al. 1996), and the Parkes-MIT- 
NRAO (PMN) radio continuum survey at 4.85 GHz (Griffith & 
Wright, 1993), as well as for sources with 6 < 0°. The FIRST 
survey used the Very Large Array (VLA) to observe the sky at 
20 cm (1.4 GHz) with a beam size of 5.4 arcsec. FIRST was 
designed to cover the same region of the sky as the SDSS, and 
observed 9000 deg 2 at the north Galactic cap and a smaller 2° .5 
wide strip along the celestial equator. It is 95 per cent complete 
to 2 mJy and 80 per cent complete to the survey limit of 1 mJy. 
The survey contains over 800,000 unique sources, with an astro- 
metric uncertainty of < 1 arcsec. 

The GB6 survey at 4.85 GHz was executed with the 91-m 
Green Bank telescope in 1986 November and 1987 October. 
Data from both epochs were assembled into a survey covering 
the 0° < 5 < 75° sky down to a limiting flux of 18 mJy, with 3.5 
arcmin resolution. GB6 contains over 75,000 sources, and has a 
positional uncertainty of about 10 arcsec at the bright end and 
about 50 arcsec for faint sources (Kimball & Ivezic 2008). The 
PMN surveys were made using the Parkes 64-m radio telescope 
at a frequency of 4850 MHz with the NRAO multibeam receiver 
mounted at the prime focus (Griffith & Wright 1993). The sur- 
veys had a spatial resolution of approximately 4'. 2 FWHM and 
were made for the southern sky between declinations of -87° 
and +10°, and all right ascensions during June and November 
in 1990. The positional accuracy is close to 10 arcsec in each 
coordinate. The survey was divided into four declination bands. 
One of these four is the equatorial survey (-9°.5 < 8 < +10°.0) 
covering 1.90 sr, which contains 11,11 A sources to a flux limit of 
40 mJy and largely overlaps the GB6 survey in the declination 
range from 0° to +10° (Griffith et al. 1995). 

The initial quasar sample was first cross-correlated between 
the SDSS quasar positions and the FIRST catalogue to within 2 
arcsec (see e.g. Ivezic et al. 2002; Lu et al. 2007). The result- 
ing sample of SDSS quasar positions was then cross-correlated 
with both the GB6 and PMN equatorial catalogues to within 1 ar- 
cmin (e.g. Kimball & Ivezic 2008). Owing to the different spatial 
resolutions of FIRST, GB6, and PMN, multiple FIRST counter- 
parts were found to within 1 arcmin for some quasars, although 
there is only single GB6 and/or PMN counterpart existed. The 
optical variability of quasars with single FIRST counterparts to 
within 1 arcmin of the SDSS positions has been presented in Gu 
& Ai (201 1). In this paper, we focus on the optical variability of 
quasars with multiple FIRST counterparts. 

There are total 29 sources having multi FIRST counterparts 
within 1 arcmin of SDSS positions. The radio spectral index a r 
was then calculated between the integrated FIRST 1 .4 GHz from 
multi-counterparts within 1 arcmin and either or both of the GB6 
and PMN 4.85 GHz. Eleven sources are defined as SSRQs with 
a x > 0.5 using the integrated 1 .4 GHz flux density. In compar- 
ison, only two sources are defined as SSRQs if only using 1.4 
GHz flux density of the closest counterpart. This implies that 
most genuine SSRQs may appear as multi-FIRST components, 
but be identified as FSRQs if using the flux density of the clos- 
est FIRST component, which is actually much smaller than that 
of radio lobes. When counterparts were found in both the GB6 
and PMN surveys, the spectral indices were consistent with each 
other. To further check the spectral index, we searched for the 
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counterparts in NRAO VLA Sky Survey (NVSS) within 1 ar- 
cmin of SDSS positions. The NVSS was also carried out using 
the VLA at 1 .4 GHz to survey the entire sky north of 6 = -40° 
and contains over 1.8 million unique detections brighter than 2.5 
mJy, however with lower spatial resolution 45 "beam -1 . Due to 
the lower spatial resolution, all sources have single NVSS coun- 
terpart, except for five sources, which have two counterparts, 
i.e. SDSS J015 105.80-003426.4, SDSS J213513. 10-005243.8, 
SDSS J221409. 96+005227.0, SDSS J220719. 77+004157. 3, and 
SDSS J233624.04+000246.0. We then re-calculated the spec- 
tral index using the flux density from single or integrated NVSS 
counterparts. We found that 25 out of 29 sources have consistent 
spectral index between using FIRST and NVSS. Three sources 
(SDSS J005905. 5 1+00065 1.6, SDSS J013352.66+01 1345.1 
and SDSS J021225.56+010056.1) changed from FSRQs to 
SSRQs using NVSS flux density, for which the part of emission 
can be likely missed in FIRST images. One souce, i.e. SDSS 
J015509.00+01 1522.5, changes from steep to flat spectral in- 
dex, but the spectral index are all close to 0.5. Therefore, it is 
still treated as SSRQ. Among the 29 sources, 14 are defined 
as SSRQs with a r > 0.5, whereas 15 objects are FSRQs with 
a r < 0.5 (see TableQ). 



2.3. Sample 

For completeness, in the following analysis we include five 
SSRQs of Gu & Ai (2011), and the resulting sample consists 
of 19 SSRQs. This sample of 19 SSRQs is listed in Table [Q 
which provides the source redshift, radio loudness, radio spec- 
tral index, black hole mass, disc bolometric luminosity, and the 
Eddington ratio. The same parameters relative to the 15 FSRQs 
are listed as a comparison. The distribution of redshift, radio 
loudness, black hole mass, and the Eddington ratio are shown in 
Fig.Q] We found that SSRQs are more extended in the distribu- 
tion of each parameters than FSRQs. The source redshift is taken 
directly from the SDSS DR7 quasar catalogue, which covers 
0.17 < z < 2.63 for SSRQs, while 0.34 <z< 2.00 for FSRQs. 
The radio loudness is from Shen et al. (201 1), which ranges from 
log R = 1.80 (SDSS J013352.66+011345.1) to log R = 5.00 
(SDSS J000622.60-000424.4) for SSRQs, whereas it is 1.55 - 
3.54 for FSRQs. However, the radio loudness was calculated 
as R = /6cm//2500, where / 6cm and / 25 oo are the flux density at 
rest-frame 6 cm and 2500A, respectively (see Shen et al. 2011 
for more details). Among 15 FSRQs, 4 sources have an inverted 
spectral index between 1.4 and 4.85 GHz with a t < 0. 

Black hole masses are estimated from the various empirical 
relations in the literature by using the luminosity and FWHM of 
broad H/3, Mg II, and C IV lines, i.e., Vestergaard & Peterson 
(2006) for H/?, and Kong et al. (2006) for Mg II and C IV. The 
luminosity and FWHM of broad Wf3, Mg II, and C IV lines 
are adopted from the measurements in Shen et al. (2011). The 
BLR luminosity Lblr is derived following Celotti, Padovani & 
Ghisellini (1997) by scaling the strong broad emission lines H/3, 
Mg II, and C IV to the quasar template spectrum of Francis et 
al. (1991), in which Lya is used as a flux reference of 100. By 
adding the contribution of Ha with a value of 77, the total rel- 
ative BLR flux is 555.77, which consists of HyS at 22, Mg II at 
34 and C IV at 63 (Francis et al. 1991; Celotti et al. 1997). From 
the BLR luminosity, we estimate the disc bolometric luminosity 
as LboI = IOLblr (Netzer 1990). It can be seen from Tableland 
Fig-Ethat the black hole mass of SSRQs ranges from 1O 8 29 M 
to 10 1008 M o for SSRQs, while 1O 844 M to 1O 9 56 M for FSRQs. 



The Eddington ratio of SSRQs log Leoi/i-Edd ranges from -2.35 
to 0.52, while in FSRQs it ranges from -2.03 to 0.30. 

3. Results 

The SDSS DR7 CAS contains the Stripe82 database, containing 
all imaging from SDSS stripe 82 along the celestial equator at 
the southern Galactic cap. It includes a total of 303 runs, cover- 
ing any given piece of the close to 270 deg 2 area approximately 
80 times. Only about one-quarter of the stripe 82 scans were ob- 
tained in photometric conditions, the remainder being taken un- 
der variable clouds and often poorer than normal seeing. For the 
runs that were non-photometric, an approximate calibration, us- 
ing the photometric frames as reference, was derived and made 
available in the CAS Stripe82 database. In this work, we directly 
use the point-spread-function magnitudes in the CAS Stripe82 
database from the photometric data obtained during the SDSS-I 
phase from data release 7 (DR7; Abazajian et al. 2009) and the 
SN survey during 2005 - 2007. The typical measurement error 
in magnitude is about 0.03 mag. 

Among 19 SSRQs, we select the sources classified as point 
sources in all observational runs. Only data with good measure- 
ments (high-quality photometry) are selected following the rec- 
ommendations in the SDSS instruction^ Moreover, we insist 
on the ugriz magnitude being brighter than the magnitude limit 
in each band, i.e. 22.0, 22.2, 22.2, 21.3, 20.5 at u, g, r, i, z, 
respectively. The data taken at cloudy conditions are also ex- 
cluded. We calculate the spectral index a v from the linear fit to 
the log f v - log v relation after applying an extinction correction 
to the ugriz flux density and taking the flux error into account. In 
most cases, the linear fit gave good fits. Each cycle of ugriz pho- 
tometry was usually completed within five minutes, i.e. quasi- 
simultaneously, therefore, the spectral index calculation will not 
be seriously influenced by any source variations. 

3.1. SSRQs 

3.1.1. Variability 

Among 19 SSRQs, SDSS J013352.66+01 1345.1 was excluded 
from our analysis owing to the low quality of the photometric 
data. All remaining 18 SSRQs show large amplitude variations 
with overall variations in r band Ar = 0.22 - 0.92 mag (see 
Table[TJ. In general, the variations in different bands show simi- 
lar trends. 

The correlation between the spectral index a v and PSF r 
magnitude was checked for all sources using the Spearman rank 
correlation analysis method. We found that 8 of 18 SSRQs show 
a significant positive correlation at a confidence level of > 99%. 
The positive correlation indicates that the source spectrum be- 
comes flatter when the source is brighter. The result shows that 
a bluer-when-brighter trend is clearly detected in these eight 
SSRQs. 

3.1.2. Broad Absorption Line Quasars 

Among the 18 SSRQs, we found that two sources are defined 
as broad absorption line quasar (BALQ) in the literature. SDSS 
J024534.07+010813.7 (z = 1.5363) was listed as a BALQ in 
Gibson et al. (2009) selected from SDSS DR5. This source is 
also included in the First Bright Quasar Survey with name of 
FBQS J0245+0108 (Becker et al. 2001). It appears as a single 

1 http://www.sdss.org/dr7/products/catalogs/flags.html 
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counterpart in the NVSS within 1 arcmin of SDSS position with 
a flux density of /nvss = 337.9 mJy and position angle of 55.5 
degree. However, it has three counterparts in the FIRST within 1 
arcmin of SDSS positions (see Fig. |2j. Three counterparts align 
with a similar position angle with that of NVSS counterpart. This 
difference is simply caused by the different spatial resolution of 
FIRST and NVSS in the case that the radio structure can be re- 
solved in former, however not in the latter. The flux density of the 
closest FIRST counterpart is 1 1 .67 mJy, while the integrated flux 
density of all counterparts is /first = 317.2 mJy. Therefore, the 
radio emission of this source is highly dominated by the large 
scale double radio lobes. The spectral index calculated for the 
integrated FIRST 1.4 GHz and PMN 4.85 GHz flux density of 
/pmn = 115 mJy is a r = 1.05, while it is a r = 0.82 when the 
GB6 4.85 GHz flux density of /gb6 = 85.7 mJy is used. In con- 
trast, the spectral index will be highly inverted if only the flux 
density of the closest FIRST counterpart is used. 

The light curves of SDSS J024534.07 +010813.7 in u, g, r, 
i, and z bands are shown in Fig. [3] The variations exhibit simi- 
lar trends in all bands, over seven years from 2000 to 2007. The 
overall variation in r band is 0.29 mag. In the observational ses- 
sions, the source became gradually fainter from 2000 to 2004, 
and then stayed quite stable from 2004 to 2007. 

In SDSS J024534.07+010813.7, a significant positive corre- 
lation present between the spectral index a v and r magnitude, 
which is shown in Fig. [4] The Spearman correlation analysis 
shows a significant positive correlation with a correlation coef- 
ficient of r s = 0.435 at confidence level of > 99.9%. The pos- 
itive correlation shows that the source spectrum becomes flat- 
ter when the source is brighter, i.e. that there is a bluer-when- 
brighter trend. With the source redshift z = 1.5363, SDSS ugriz 
wavebands correspond to the wavelength range of 1400 - 3521 A 
in the source rest-frame. For a sample of 17 radio-quiet AGNs, 
Shang et al. (2005) show that the spectral break occurs at around 
1 100 A for most objects, which is thought to be closely related 
to the big blue bump. If this spectral break also exists in SDSS 
J024534.07+010813.7, we would expect to observe the rising 
part of accretion disk thermal emission when it dominates over 
the nonthermal emission. However, the spectral index in the ob- 
servational sessions are in the range of 1.59 to 2.03. This implies 
a falling SEDs in the rest frame 1400-3521A. This can be likely 
caused by the internal extinction in quasar itself, for example, 
from the outflowing absorption gas. 

We found another BALQ SDSS J021728.62-005227.2 (z = 
2.4621) in SSRQs of Gu & Ai (2011). It is defined as BALQ 
in Trump et al. (2006) selected from SDSS DR3. It is defined 
as SSRQs in Gu & Ai (201 1) with radio spectral index of a r = 
0.75. Its overall variations amplitude in r band is Ar = 0.35. 
The significant positive correlation between the r magnitude and 
the optical spectral index shows a bluer-when-brighter trend in 
this source. However, unlike SDSS J024534.07+010813.7, the 
spectral index varies from 0.45 to 1.28, with most time staying 
at < 1 .0. The internal extinction from the outflowing absorption 
gas may not be severe in this source. 

3.2. FSRQs 

All 15 FSRQs show large amplitude variations with overall vari- 
ations in r band Ar = 0.18 - 0.97 (see Table Q]). The correla- 
tion between the spectral index a v and PSF r magnitude was 
checked for all sources using the Spearman rank correlation 
analysis method. We found that 11 of 15 FSRQs show a sig- 



nificant positive correlation at a confidence level of > 99%, i.e. 
a bluer- when-brighter trend. 

3.3. BLR geometry 

Assuming the radio spectral index as an indicator of orientation, 
the BLR geometry can be investigated when line width of broad 
emission lines are included (Fine et al. 2011). To this end, we 
combine all sources in this work and those FSRQs in Gu & 
Ai (2011), resulting a sample of 62 quasars, which includes 18 
SSRQs and 44 FSRQs. The relationship between the radio spec- 
tral index and FWHM of Mg II line is shown in Fig. [5] for a 
subsample of 55 sources (42 FSRQs and 13 SSRQs) with avail- 
able Mg II line measurements. There is a trend of broader Mg II 
lines with steeper spectral index. The Spearman rank correlation 
analysis shows a positive correlation with correlation coefficient 
of r s = 0.298 at confidence level of 97.3%. The mean value of 
Mg II FWHM is 4075 km s 1 for FSRQs, while it is 61 10 km s 1 
for SSRQs. The median value of Mg II FWHM is 3586 km s _1 
and 5694 km s _1 for FSRQs and SSRQs, respectively. Although 
our sample is much smaller than theirs, the results are consis- 
tent with that of Fine et al. (201 1), which used the radio spectral 
index between Westerbork Northern Sky Survey 330 MHz and 
NVSS 1.4 GHz. The result implies a disk-like BLR geometry 
may present in these quasars. 

4. Discussions 

4.1. SSRQs and FSRQs 

While FSRQs are usually associated with core-dominated radio 
quasars, SSRQs are generally related to lobe-dominated ones, 
usually with two large-scale optically thin radio lobes. Usually, 
the beaming effect is not severe in SSRQs because of the rela- 
tively large viewing angle. Therefore, the optical continuum of 
SSRQs could be dominated by the thermal emission. In Fig. [6] 
most SSRQs with available Mg II measurements lie close to or 
below the solid line for radio quiet AGNs. The average value of 
the ratio of the luminosity at 3000 A to that of broad Mg II line 
is < log (vL v (3000A)/L M gii) >= 1-78 + 0.26, which is in good 
agreement with the linear relation for radio quiet AGNs (see Fig. 
|6}. This likely implies that the thermal emission from accretion 
disk is indeed the dominant one in the optical continuum at the 
epoch of SDSS spectra for most SSRQs, if not all. For compari- 
son, 14 FSRQs with Mg II measurements are also plotted in Fig. 
[6] We found that the distribution of FSRQs are generally similar 
to that of SSRQs, but with larger scatter. The average value of 
the ratio of the luminosity at 3000 A to that of broad Mg II line 
is < log (vL v (3000A)/L MgII ) >= 1.81 + 0.31, which is slightly 
larger than that of SSRQs. This may imply that the optical non- 
thermal jet emission may not be dominant in all FSRQs, as found 
by Chen et al. (2009). 

Interestingly, eight of 18 SSRQs display a BWB trend (see 
Table [TJ. The BWB trend has also been found in radio quiet 
AGNs, and the variability was found to be anti-correlated with 
the Eddington ratio (Ai et al. 2010). In radio quiet AGNs, the 
optical variability could be due to the variation in the accre- 
tion process, for example, the variation in the accretion rate 
(e.g. Li & Cao 2008), since the jet is either weak or absent. 
These explanations could also be used to explain the variabil- 
ity of SSRQs if the optical emission is indeed mainly from the 
accretion disk. To further investigate the mechanism of variabil- 
ity in SSRQs, the relationship between the variability Ar and the 
Eddington ratio is plotted in Fig. [7] We found an anti-correlation 
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with a Spearman correlation coefficient of r s = -0.490 at confi- 
dence level of 96.2%. This similar anti-correlation to radio quiet 
AGNs, strongly implies that the mechanisms of optical variabil- 
ity may be similar in the two populations. The thermal emis- 
sion from the accretion disk may be responsible for the vari- 
ability of SSRQs. In contrast, there is no correlation for FSRQs. 
This implies that the mechanism of variability in FSRQs may 
be different from that of SSRQs. However, that the variation in 
SSRQs is caused by jet nonthermal emission cannot be com- 
pletely excluded in some cases because the variability ampli- 
tudes are generally larger that the typical values for radio quiet 
AGNs, 0.05 - 0.3 mag (e.g. Ai et al. 2010), (see TableQ]). In ex- 
treme cases, the overall amplitude is Ar = 0.92 mag in SDSS 
J013514.39-000703.8. It may be more likely that both thermal 
and nonthermal emission contribute to the variability. Further 
multi-waveband monitoring is needed to help resolve these un- 
certainties, especially the spectroscopic monitoring. 

In terms of optical variability on month-to-year time-scales, 
no obvious difference was found in the behavior of radio quiet 
quasars and lobe-dominated quasars (LDQs) (Stalin et al. 2004). 
Moreover, there is no significant difference in either the ampli- 
tude or duty cycle of intranight optical variability (INOV) be- 
tween these two classes of non-blazar AGNs. However, their 
sample is too small, and the time spanned by observations is 
too short. Nevertheless, the authors infer that the radio loudness 
of a quasar alone is not a sufficient condition for a pronounced 
INOV. While the observational data cannot exclude accretion 
disc flares as the source of the intranight optical variability in 
radio quiet quasars and LDQs, it does not preclude a substan- 
tial contribution from blazar-like relativistically beamed emis- 
sion. As suggested by Stalin et al. (2005), the radio quiet quasars 
may also eject relativistic jets. However, their jets are probably 
quenched while crossing the innermost micro-arcsecond scale, 
possibly through heavy inverse Compton losses in the vicinity 
of the central engine (e.g. Gopal-Krishna et al. 2003). However, 
the analysis of spectra of a set of radio quiet quasars that had 
already been searched for microvariability showed that the jet- 
based scenario is unlikely (Chand, Wiita & Gupta 2010). 

In order to compare systematically the variations of SSRQs 
with FSRQs, we combine the FSRQs of Gu & Ai (2011) with 
sources in this work. The combined sample consists of 62 
quasars, of which 44 are FSRQs, and 18 are SSRQs. We show 
the histogram of variation amplitude in r band Ar for both 
SSRQs and FSRQs in Fig. [8] It can be seen that Ar of FSRQs 
are systematically larger than that of SSRQs, with median value 
of 0.52 mag and 0.41 mag for FSRQs and SSRQs, respectively. 
Apparently, FSRQs cover wider range in Ar. Four FSRQs ex- 
hibit Ar > 1 .0 mag with an extreme case of Ar = 3.46 (see Gu & 
Ai 201 1), while none in SSRQs. This implies that the variations 
in FSRQs are systematically more violent than SSRQs. This is 
most likely due to the relativistic jet, which is viewed at smaller 
viewing angle in FSRQs. The fact that the typical variations am- 
plitude of SSRQs is between the values of radio quiet AGNs 
and FSRQs is qualitatively consistent with the source nature of 
SSRQs, which are radio-loud sources with prominent powerful 
relativistic jets, however viewed at larger viewing angle resulting 
in a dominated thermal emission at optical wavebands. In terms 
of the spectral variations, eight of 18 SSRQs show BWB trend 
at confidence level of > 99%, i.e. at probability of ~ 44%. In 
contrast, 25 (~ 57%) of the 44 FSRQs show BWB trend at con- 
fidence level of > 99%, while only one source shows a RWB 
trend at confidence level of > 99%. Therefore, it seems that 
BWB trend are quite common both in FSRQs and in SSRQs, al- 
though the mechanisms can be different, for example, shock-in- 



jet and accretion rate change are mainly responsible for FSRQs 
and SSRQs, respectively. 

In the unification scheme of AGNs, the radio galaxies are 
generally unified with radio quasars (Antonucci 1993), with ra- 
dio galaxies being viewed at large viewing angle. The jet emis- 
sion in radio galaxies are not severely boosted, then the dom- 
inated optical emission is expected to be thermal emission re- 
lated to the accretion disk. Indeed, it is argued that the optical 
and X-ray emission arise from the accretion disk - corona sys- 
tem in two radio galaxies 3C 120 and 3C 1 1 1 by using extensive 
multi-frequency monitoring data (Chatterjee et al. 2009, 2011). 
In both sources, the low optical polarization supports the thermal 
origin of the optical emission. However, in 3C 111 the detection 
of polarization percentages of 3% or even higher suggests that 
during particular activity phases an additional contribution from 
non-thermal emission from the jet may be present (Jorstad et al. 
2007; Chatterjee et al. 201 1). Interestingly, signicant dips in the 
X-ray light curve are followed by ejections of bright superlu- 
minal knots in the VLBA images for both sources. The authors 
argued that the radiative state of accretion disk plus corona sys- 
tem, where the X-rays are produced, has a direct effect on the 
events in the jet, where the radio emission originates. 

Similar to many FSRQs in Gu & Ai (201 1) having a spectral 
index a v < 1 .0, almost all FSRQs in this work show a v < 1 .0, 
implying a rising SED in optical regions. One possibility is that 
the optical emission is mainly from thermal accretion disc, al- 
though these sources are defined as FSRQs. Alternatively, they 
can be high-synchrotron-peaked FSRQs, since the synchrotron 
peak frequency is higher than the rest-frame frequency of SDSS 
ugriz wavebands. This however needs further investigations. 

4.2. BALQs 

The nature of BAL quasars is still unclear. It was argued that 
the outflowing BAL wind is preferentially equatorial, and only 
in those objects with almost edge-on accretion disk to the line 
of sight, can be observed as BAL (Cohen et al. 1995; Goodrich 
& Miller 1995). However, the VLBI radio images of BALQs 
argued against this orientation scenario, in terms of the radio 
structures and both the steep and flat radio spectrum (e.g Jiang 
& Wang 2003; Liu et al. 2008; Doi et al. 2009). Fine et al. 
(201 1) found no evidence that BAL QSOs have a different spec- 
tral index distribution to non-BALs although only 25 obvious 
BALs are considered in their sample. Their BALQs both have 
steep and flat radio spectrum computed using the Westerbork 
Northern Sky Survey (WENSS) at 330 MHz and the NVSS at 
1.4 GHz. In general, the radio-loud BALQs tend to be com- 
pact in the radio, similar to Gigahertz-Peaked Spectrum (GPS) 
and Compact Steep Spectrum (CSS) sources, which are thought 
to be at the young stage of large scale radio sources (Becker 
et al. 2000; Montenegro-Montes et al. 2009). It has been pro- 
posed that the BALs are not closely related with the inclina- 
tions, and may be associated with a relatively short-lived evo- 
lutionary phase with a large BAL wind covering fraction (e.g. 
Briggs et al. 1984; Gregg et al. 2000). We found that the two 
radio-loud steep spectrum BALQs in this paper have different 
radio structures. SDSS J024534.07+010813.7has prominent ex- 
tended double radio lobes, of which the emission dominates 
over that of radio core (see Fig. [2]). The source is too extended 
(~ 52 arcsec, i.e. ~ 440 kpc) to be considered as a CSS/GPS. 
This source has not been imaged with VLBI. In contrast, SDSS 
J021728.62-005227.2 is compact in FIRST image with only 
single counterpart found. It has been imaged with VLBI by Doi 
et al. (2009). However, only a compact component was found 
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with flux density at 8.4 GHz /g.4GHz = 43 mJy at scale of < 1 
kpc. In combination with the FIRST flux density, the radio spec- 
trum is steep a r - 0.9, which makes this source resemble to CSS, 
thus likely being a young radio source. More radio observations 
at higher spatial resolution would be helpful to better understand 
the nature of this source. 

While there are many investigations on the variability of 
broad absorption line troughs (e.g. Barlow 1993; Lundgren et 
al. 2007; Gibson et al. 2008, 2010; Capellupo et al. 2011), not 
much works have been done on the optical continuum variabil- 
ity for BALQs. On timescales of the order of years, the optical 
continuum variability of BALQs was found to differ markedly 
from that of non-BAL radio quiet quasars (Turnshek 1988). In 
short timescales, intranight optical variations of ~ 5% - 9% 
were detected on timescales of ~ 1 hr in two BALQs (Anupama 
& Chokshi 1998). The recent search for optical microvariabil- 
ity in a sample of 19 radio quiet BALQs argued that radio quiet 
BALQs do not appear to be a special case of the radio quiet 
quasars in terms of their similar duty cycle (Joshi et al. 2011). 
However, all these observations focused only on the radio quiet 
BALQs. Therefore, for the first time, we present here the opti- 
cal continuum variability for radio loud BALQs, although only 
two sources are found in our sample. In both BALQs, the C IV 
BAL troughs fall outside the ugriz bands. Together with the sim- 
ilar trends of variability in different wavebands rule out the pos- 
sibility of the observed light-curve variations being caused by 
variations in the strength of the BAL features. The variation am- 
plitude in the r band of two sources is not atypical in overall 
sample of SSRQs, although they are at the low end of distribu- 
tion in Fig. [8] This implies that the mechanism of optical con- 
tinuum variability of radio-loud BALQs may likely be similar 
to that of more general SSRQs. Interestingly, both BALQs show 
BWB trend. The percentage of BWB trend is higher than that of 
SSRQs. However, the small sample size precludes to draw any 
firm conclusions. Certainly, a larger sample is needed for further 
investigations. The spectral index of two BALQs in this work 
are steeper than the typical value of other SSRQs, which implies 
more reddened in BALQs than non-BALQs. This is consistent 
with previous results (Trump et al. 2006; Gibson et al. 2009). 

5. Summary 

We have constructed a sample of 19 SSRQs and 15 FSRQs 
in the SDSS Stripe 82 region. The variability and the relation- 
ship between the spectral index and brightness were investigated 
for 18 SSRQs and 15 FSRQs. We found that all SSRQs show 
large-amplitude overall variations, e.g. from 0.22 to 0.92 mag 
in r band. We found a bluer-when-brighter trend in 8 of the 18 
SSRQs studied here. In contrast, the variability amplitude of 15 
FSRQs in r band ranges from 0.18 to 0.97 mag, and 11 of 15 
sources show a BWB trend. We found an anti-correlation be- 
tween the Eddington ratio and the variability amplitude in r band 
for SSRQs, which is similar to that in radio quiet AGNs. This 
implies that the thermal emission from accretion disk may be re- 
sponsible for the variability in SSRQs, although the jet nonther- 
mal emission cannot be excluded. In contrast, there is no cor- 
relation in FSRQs, which implies that the mechanisms of vari- 
ability in FSRQs may be different from that in SSRQs. The opti- 
cal continuum variability of radio loud BALQs is investigated 
for the first time here on two sources with steep radio spec- 
trum. Both radio loud BALQs show variations with amplitude 
of about 0.3 mag at r band. Their spectral variability all show a 
bluer-when-brighter trend. In combination with 29 FSRQs also 
selected from Stripe 82 region in Gu & Ai (2011), we found a 



trend of a broader line width of broad Mg II line with steeper ra- 
dio spectral index. This implies a disc-like BLR geometry may 
present in these quasars. In these 62 quasars, we found a ~ 57% 
source percentage showing BWB trend in FSRQs, whereas it is 
~ 44% in SSRQs. 
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Fig. 1. The histogram of sample parameters: redshift (upper left), 
radio loudness (upper right), black hole mass (lower left), and 
the Eddington ratio (lower right). The dotted lines are for 19 
SSRQs, and solid lines are for 15 FSRQs. 
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Fig. 2. The FIRST radio image of SDSS J024534.07+010813.7. 
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Fig. 3. The ugriz band light curve 
J024534.07+010813.7 (from top to bottom). 
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Fig. 4. The relationship between the spectral index and the PSF 
magnitude at r band for SDSS J024534.07 +010813.7. A signifi- 
cant positive correlation is present, which implies a bluer-when- 
brighter trend. 
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Fig. 5. The line width of broad Mg II line versus radio spectral Fig. 8. The histogram of Ar for SSRQs and FSRQs in Gu & Ai 
index. The solid circles represent SSRQs, while the open circles (201 1) and in this work. The dotted lines are for 18 SSRQs, and 
are for FSRQs. solid lines are for 44 FSRQs. 
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Fig. 6. The plot of broad Mg II line luminosity and continuum 
luminosity at 3000 A for 13 SSRQs and 14 FSRQs. The sym- 
bols are same as in Fig. 4. The square is for BALQ SDSS 
J024534.07+010813.7. The solid line is the OLS bisector lin- 
ear fit to radio-quiet AGNs in Kong et al. (2006), AL A 3000 ^ = 
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Fig. 7. The Eddington ratio versus the variability at r band Ar for 
SSRQs. 
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Table 1. Source list: Col. 1 - SDSS source name; Col. 2 - red- 
shift; Col. 3 - radio loudness; Col. 4 - the spectral index between 
1 .4 and 4.85 GHz, of which GB6 flux density is used for those 
labeled with a , otherwise PMN one is used; Col. 5 - black hole 
mass in unit of solar mass; Col. 6 - disc bolometric luminosity in 
unit of erg s _1 ; Col. 7 - the Eddington ratio / = Lbol/^edd; Col. 
8 - overall variation in r band; Col. 9 - the range of variation in 
the spectral index a v ; Cols. (10-11) - the Spearman correlation 
coefficient and probability level, respectively. 
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from 


Gu&Ai(2011) 












J012401.76+003500.9 


1.8516 


3.85 


0.86" 


9.36 


46.11 


-1.36 


0.84 


-0.52,0.61 


0.687 


6.8e-07 


J012517.14-001828.9 


2.2780 


3.44 


0.59 


8.50 


47.13 


0.52 


0.29 


0.22,0.87 


-0.202 


0.11 


J015832.51-004238.2 


2.6071 


3.58 


0.87 


8.29 


46.90 


0.50 


0.47 


0.44,1.19 


-0.207 


0.28 


J021728.62-005227.2 


2.4621 


3.29 


0.75 


8.84 


46.40 


-0.55 


0.35 


0.45,1.28 


0.488 


7.7e-06 


J022508.07+001707.2 


0.5270 


3.81 


0.86" 


8.88 


45.19 


-1.80 


0.53 


0.25,0.68 


0.470 


0.03 


FSRQs 


J000111. 19-002011.5 


0.5179 


2.58 


0.15 


8.57 


45.16 


-1.52 


0.36 


0.64,1.04 


0.426 


5.44e-03 


J010033.50+002200.1 


0.7534 


1.94 


0.11" 


9.21 


45.93 


-1.40 


0.55 


-0.33,0.26 


0.615 


2.39e-08 


J010826.84-003724.2 


1.3724 


3.31 


0.26 


8.76 


46.31 


-0.56 


0.65 


0.25,0.80 


0.330 


4.00e-03 


J015105.80-003426.4 


0.3352 


1.55 


-0.24 


8.73 


44.82 


-2.03 


0.55 


0.05,0.95 


0.807 


2.39e-12 


J020234.32+000301.7 


0.3664 


3.26 


0.31 


8.47 


44.81 


-1.77 


0.18 


0.19,0.58 


-0.274 


0.36 


J020326.98+003744.3 


1.5840 


3.43 


0.30" 


8.44 


46.31 


-0.24 


0.80 


0.24,0.89 


0.678 


2.71e-05 


J021612.20-010518.9 


1.4931 


2.61 


0.03 


8.77 


46.45 


-0.43 


0.27 


0.37,0.73 


0.263 


0.04 


J024854.81+001053.8 


1.1457 


2.92 


0.30 


8.65 


46.00 


-0.77 


0.97 


-0.31,0.56 


0.723 


7.06e-12 


J025928.51-001959.9 


2.0001 


2.54 


-0.63 


8.87 


47.29 


0.30 


0.30 


0.08,0.66 


0.323 


1.24e-02 


J031318.66+003623.9 


1.2561 


2.35 


0.11" 


9.56 


46.02 


-1.65 


0.32 


0.52,0.98 


0.353 


4.51e-03 


J220719.77+004157.3 


1.8926 


3.36 


-0.05 


8.85 


45.90 


-1.06 


0.59 


0.05,0.58 


-0.171 


0.18 


J222729.05+000521.9 


1.5133 


3.02 


0.15 


9.19 


46.52 


-0.79 


0.70 


0.46,1.36 


0.722 


3.18e-10 


J233200.00+01 1510.9 


1.6395 


3.45 


0.27" 


8.75 


46.01 


-0.86 


0.44 


0.15,0.87 


0.647 


8.97e-08 


J234624.56+001914.2 


1.7775 


2.84 


-0.02" 


8.87 


46.39 


-0.60 


0.79 


0.09,0.72 


0.702 


9.77e-ll 


J234830.98+01 1037.6 


1.7064 


3.54 


0.19" 


8.44 


45.87 


-0.68 


0.60 


0.17,1.08 


0.746 


1.67e-07 



